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MILC collaboration results for /b, /s^ , fo, /csand their ratios are presented. These results are still preliminary, 
but the analysis is close to being completed. Sources of systematic error, both within the quenched approximation 
and from quenching itself, are estimated. We find, for example, fs ~ 153 ±10 '^q^ MeV, and fB^/fs = 1.10 ± 
0.02 lo o3 lo 021 where the errors are statistical, systematic (within the quenched approximation), and systematic 
(of quenching), respectively. The extrapolation to the continuum and the chiral extrapolation are the largest 
sources of error. Present central values are based on linear chiral extrapolations; a shift to quadratic extrapolations 
would raise /s by Ri20 MeV and make the error within the quenched approximation more symmetric. 



The MILC collaboration is continuing its pro- 
gram of calculating the decay constants of 
heavy-light pseudoscalar mesons. The computa- 
tions use Wilson light quarks and Wilson and 
static heavy quarks. We work on both quenched 
lattices, with a wide range of lattice spacings, and 
Np = 2 dynamical staggered lattices. Table || 
gives the lattice parameters. 

A major improvement in the past year has 
been the completion of dedicated runs to de- 
termine the static-light decay constants on lat- 
tices A,B,Q,E,G,L,N,0,M,P. These runs use a 
multi-source technique, with relative wavefunc- 
tions taken from ||]. On lattices C,D,H,F,G 
— generally the smaller physical volumes — we 
get acceptable static-light data as a simple by- 



product of the Henty-Kenway hopping parameter 
expansion 1^ used for the heavy Wilson quarks. 
On lattice G, where both methods are available, 
the results are consistent. 

A second improvement has been the calcula- 
tion, following the approach of [Q, of the scale q* 
of the coupling in the perturbative renormaliza- 
tion constant Za of the axial current. For propa- 
gating Wilson quarks, the result is, after tadpole 
improvement, q* — 2.32/ a Mass dependent 
effects are not included at this point. We esti- 
mate the systematic error of the renormalization 
by changing q* by a factor of 2 and reanalyzing. 
The error is rather small (<3%). 

Currently, we find (in MeV): 
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/B = 153(10)(^t^)(+i^); /b. 



164(9)(itD(±o') 



The errors are statistical (plus "fitting"), system- 
atic (within the quenched approximation), and 
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Table 1 

Lattice parameters. Lattices F, G, and L-P use 
variable-mass Wilson valence quarks and two fla- 
vors of fixed-mass staggered dynamical fcrmions; 
all other runs use quenched Wilson quarks. Lat- 
tice G was generated by HEMCGG; lattice F, by 
the Columbia group. 
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A 


5.7 


8^ X 48 


200 


B 


5.7 


16^ X 48 


100 


E 


5.85 


12^ X 48 


100 


Q 


6.0 


12^ X 48 


235 


c 


6.0 


16^ X 48 


100 


D 


6.3 


24^ X 80 


100 


H 


6.52 


32^ X 100 


60 


F 


5.7 (0.01) 


16^ X 32 


49 


G 


5.6 (0.01) 


16^ X 32 


200 


L 


5.445 (0.025) 


16^ X 48 


100 


N 


5.5 (0.1) 


243 X 64 


100 





5.5 (0.05) 


243 X 64 


100 


M 


5.5 (0.025) 


20^ X 64 


100(200) 


P 


5.5 (0.0125) 


20^ X 64 


199 



systematic (of quenching), respectively. For the 
ratios, we get: 

fBjfB = 1.10(2)(t^)(+3) 

fnJfD - im{2){tl)ei) 

fs/fD, = 0.76(3)(tD(to)- 

A discussion of the most important sources of sys- 
tematic errors follows. 

• Chiral Extrapolation. The proper functional 
form to use in extrapolating physical quantities 
to the chiral limit is not clear. For example, al- 
though lowest order chiral perturbation theory 
predicts that is a linear function of quark 
mass, we observe ||^] small but significant de- 
viations from linearity. 

These deviations could be due to unphysical 
effects such as the finite lattice spacing or vol- 
ume. In addition, the curvature can be changed 
significantly, and sometimes made negligible, by 
shifting the fit range (in t) on the individual prop- 
agators. Even the more "physical" cause (chiral 
logs or higher order analytic terms) are a source 



of possible spurious effects because quenched chi- 
ral logs are in general different from those in the 
fuU theory §. 

For these reasons, we presently fit quantities 
like to their lowest order chiral form, despite 
the poor confidence levels. The systematic er- 
ror is estimated by repeating the analysis with 
quadratic (constrained) fits. This error is < 10% 
for decay constants on all quenched data sets used 
to extrapolate to the continuum; usually it is 
<5%. (After extrapolation to the continuum, the 
error is larger: 7% to 15%.) 

Our reasons for choosing linear chiral fits for 
the central values are somewhat subjective, and 
it is possible that we will switch to quadratic 
fits in the final version of this work. To help us 
make the choice, we are studying a large sample of 
quenched lattices at /3 = 5.7, with volumes up to 
24^^ On this sample we have six light quark 
masses (as opposed to three for each of the lattices 
used for the heavy-light computation) and have 
light-light mesons with nondegenerate as well as 
degenerate quarks. 

If a switch to quadratic fits were to be made 
now, it would raise the central values of fs, fs,, 
fo and fn^ by 23, 19, 13, and 14 MeV, respec- 
tively. The systematic error within the quenched 
approximation would then become much more 
symmetric, with the continuum extrapolation the 
dominant positive error and the chiral extrapola- 
tion the dominant negative one. 

• Heavy-Quark Interpolation. Having static- 
light results on all lattices allows us to find decay 
constants for physical B mesons by interpolation 
between heavy-light and static-light data, rather 
than extrapolation from the former. The interpo- 
lating fits have good confidence levels on all our 
data sets, reassuring us that the procedure pT| | 
we use for the heavy-light data is reasonable. 

One estimate of the systematic error of this ap- 
proach is obtained by comparing decay constants 
computed with two different mass ranges of prop- 
agating quarks: "lighter heavies," (mesons 1.25 
to 2 GeV) and "heavier heavies," (mesons 2 to 4 
GeV). The difference is less than 1% at the three 
weakest quenched couplings (/3 = 6.0,6.3,6.52), 
less than 5% over all lattices, and less than 4% 
after linear extrapolation of all quenched lattices 
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to a = 0. 

For quenched lattices A,B,E, the new static- 
Ught resuhs produce only small changes from that 
reported previously However, on the large 
Nf = 2 lattices, including the static point raises 
Jb by about one old (statistical) standard devia- 
tion, and reduces the statistical error (and the dif- 
ference between using heavier-heavies and lighter- 
heavies) by about 50%. 

• Extrapolation to the Continuum. For any 
physical quantity Q computed here, we expect 
Q{a) = (5a=o(l + aMi -!-•••). In practice, we 
find the slope to be large for the decay constants 
{Ml - 300-650 MeV), with /s, the worst of- 
fender. This leads to large extrapolation errors 
(~ 12-27%). The ratios of decay constant are 
much better behaved, with A4i ^ 100 MeV and 
an error of ~4-5%. 

Figure |^ shows several fits of Jb vs. a used to 
estimate the two largest sources of systematic er- 
ror. The central value is obtained from a linear fit 
to all the diamonds, which in turn use linear chi- 
ral fits, a lattice scale set by /jr, and the "EKM" 
corrections . The error of the continuum ex- 
trapolation is estimated by comparing the central 
value with the results of a constant fit to the three 
diamonds with smallest values of a, a linear fit to 
the octagons (which use "lighter heavies" and no 
EKM corrections) , and two other similar types of 
fits (not shown). The continuum extrapolation 
error is defined as the largest of these four differ- 
ences and is in practice almost always given by the 
first difference. The difference of the extrapola- 
tion of the squares (which have a quadratic chiral 
extrapolation) and the central value determines 
the chiral extrapolation error. 

• Quenching Effects. We have repeated our 
computations on several Nf—2 dynamical fer- 
mion lattices (crosses in Fig. Such compu- 
tations are not yet "full QCD" because (1) the 
virtual quark mass is fixed and not extrapolated 
to the chiral limit, (2) the Np — 2 data is not 
yet good enough to extrapolate to a — 0, and (3) 
we have two light flavors, not three. Thus the 
Np = 2 simulations are used at this point only for 
systematic error estimation. See ^ for details. 

We thank CCS (ORNL), Indiana University, 
SDSC, PSC, CTC, MHPCC, CHPC (Utah), and 
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Figure 1. Results for /b as a function of lattice 
spacing. The diamonds are quenched approxima- 
tion points. Octagons and squares are used for 
estimating systematic errors — see text. Crosses 
are from Np = 2 lattices but are otherwise com- 
puted like the diamonds. 
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